High fidelity measurement and quantum feedback control in circuit QED 
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Circuit QED is a promising solid-state quantum computing architecture. It also has excellent 
potential as a platform for quantum control - especially quantum feedback control - experiments. 
However, the current scheme for measurement in circuit QED is low efficiency and has low signal- 
to-noise ratio for single shot measurements. The low quality of this measurement makes the imple- 
mentation of feedback difficult, and here we propose two schemes for measurement in circuit QED 
architectures that can significantly improve signal-to-noise, and potentially achieve quantum limited 
measurement. Such measurements would enable the implementation of quantum feedback protocols 
and we illustrate this with a simple entanglement stabilization scheme. 
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I. INTRODUCTION 

The emergence of quantum computing has revitalized 
the field of quantum control Q, 0, y, in recent years. 
The prospects for a working quantum computer are inti- 
mately dependent on the ability to reliably control indi- 
vidual and coupled quantum systems. There have been 
several experimental proposals for architectures capable 
of quantum computing [5j and most are in their initial 
stages of experimental development. During these ini- 
tial stages, quantum computing experiments are cither 
not possible or only proof of principle - mainly due to 
the small number of qubits under control. However, 
many of these quantum computing architectures, dur- 
ing their initial stages of experimental development, are 
ideal testbeds for quantum control experiments designed 
to investigate the limits of quantum control, quantum 
chaos, and the fundamentals of quantum mechanics. 

A particularly promising group of quantum computing 
architectures that are well suited to quantum control ex- 
periments are the solid-state realizations of cavity quan- 
tum electrodynamics. These architectures couple super- 
conducting electronic circuit elements, which serve as the 
qubits, to either electromagnetic modes 0,0 or nanoelec- 
tromechanical resonators [8|, which serve as a 'quantum 
bus' that mediates inter-qubit coupling and/or facilitates 
measurement of a qubit's state. Strong coupling between 
the qubits and the quantum bus, and the substantial 
amount of control afforded by these architectures over 
both subsystems makes these systems excellent quantum 
control testbeds. 

In this paper we will focus on one such implementation, 
that of Schoelkopf et. al. Q at Yale University, and 
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examine it from a quantum control perspective. This 
implementation couples the quantized charge degree of 
freedom of a Cooper pair box (CPB) yj, UJj, |Tl| to a har- 
monic oscillator mode of a microwave resonator. Both 
the Cooper pair box and the resonator are fabricated on 
an integrated circuit and thus the setup has been dubbed 
'circuit QED'. A number of features of this architecture 
make it well suited to quantum control - especially quan- 
tum feedback control experiments: the fact that the 
CPBs (which can be viewed as artificial atoms) are local- 
ized and the high degree of electrical control over them, 
the strong coupling between the CPBs and the resonator 
which allows for a well controlled 'atom'-field interaction, 
the fact that feedback can be applied to either the CPB 
or the microwave resonator mode, and of course, the in- 
tegrated nature of the whole apparatus. However, the 
current scheme for measurement of the qubits is highly 
inefficient and has a poor signal to noise ratio (for single 
shot measurements) |l2j . The primary reason for this 
is that while the qubit and resonator are at a tempera- 
ture of 20mK, the mixing electronics and local oscillator 
needed for amplitude measurements are at room temper- 
ature. See figure 2 in Ref. @- The noise introduced 
by this room temperature mixing, and the amplification 
stage before it, results in a measurement signal with low 
single-shot signal-to-noise ratio. This low fidelity mea- 
surement scheme is a problem from a quantum control 
perspective because high quality measurements are es- 
sential to quantum control, especially quantum feedback 
control, experiments. Ideally the measurements should 
be quantum-limited for quantum feedback control. Mo- 
tivated by this, we propose two methods for high fidelity 
measurement in the circuit QED architecture. The first 
uses an single electron transistor (SET) as a heterodyne 
mixer and the second performs the local oscillator mixing 
at the cryogenic level by the addition of a second trans- 
mission line resonator that carries a local oscillator to the 
circuit . 

This paper is organized as follows: sections [H] and 
IIHI review the circuit QED architecture of Wallraff et. 
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at, and techniques for SET based measurement. Sec- 
tion ]W\ presents our two proposals for high fidelity mea- 
surements in circuit QED. Section describes a simple 
feedback protocol for deterministically creating and sta- 
bilizing entanglement between two qubits that would be 
made possible by such high fidelity measurements, and 
finally, section IVll concludes with a discussion. 




II. CIRCUIT QED 



The qubits in the circuit QED architecture are split 
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junction Cooper pair boxes These devices can 

be modelled as two-level systems with the Hamiltonian 



(1) 



where E e i is the electrostatic energy and Ej is the 
Josephson coupling energy. The Pauli matrices d z and 
o x are in the charge basis; that is, the basis states cor- 
respond to either zero or one excess Cooper-pair charges 
on the island. 

These CPBs can be viewed as artificial atoms with 
large dipole moments, and in circuit QED they are cou- 
pled to microwave frequency photons in a quasi-one di- 
mensional transmission line cavity (a co-planar waveg- 
uide resonator) by an electric dipole interaction 6] . This 
apparatus has a number of in situ tunable parameters, 
including E e i and Ej, and a choice can be made 01 
such that the combined Hamiltonian, for qubit and trans- 
mission line cavity, is the well-known Jaynes-Cummings 
Hamiltonian flFfl: 
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H = haj r (a'a H — ) H a z + hg(a'a 



(2) 



where a (a^) is the annihilation (creation) operator for 
the cavity mode, uj t is the cavity resonance frequency, ui a 
is the energy splitting of the qubit, and g is the coupling 
strength. Of course, the qubit energy splitting, u> a , is 
a function of the Cooper pair box parameters E e i and 
Ej 0|- For typical values of the parameters in this 
Hamiltonian see Refs. |(| It should be noted here 
that there has been a basis change between Eqns. 
and In essence, we have swapped the x and z axes 
of Eq. and the computational basis for our qubit has 
become the Josephson basis of the CPB. This choice of 
basis, which corresponds to operating the CPB in what 
is called the charge degeneracy point {E e i ~ 0), has a 
number of advantages, the primary one being that the 
computational basis states become first-order insensitive 
to dephasing from offset charge noise In fact, a CPB 
is only effective as a robust qubit at this operating point. 

In the dispersive regime, where A = uj a — oj r ^> 1, the 
above Hamiltonian can be diagonalized to obtain [lJI 



H m h(u) r + 
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FIG. 1: Modeling the resonator as a two-sided leaky cavity. 



This form of the Hamiltonian makes evident the qubit- 
state dependent pull of the cavity frequency in this 
regime. Thus readout of the qubits in this architecture 
can be done by irradiating the cavity and then probing 
the reflected or transmitted amplitude. A quantum non- 
demolition (QND) measurement of the qubit in the a z 
basis can be made by measuring the phase shift (f> of a 
probe microwave through the cavity at frequency u) r |l4| . 
Physically, this measurement localizes the qubit in the 
Josephson basis, the eigenbasis of <j z in Eq. @, and the 
CPB remains operating at the optimal charge degeneracy 
point. 

The precise dependence of the phase of the output field 
on the qubit state can be seen from the input-output 
relations between the in and out cavity fields. If we model 
the resonator as a two-sided cavity with equal rates of 
leakage on both ends, as in figure^ then the relationship 
between the mode in the cavity and the output mode is 

M 



a ou t(t) = y/Ka(t) 



(4) 



where a ou t is the output mode, a is the intra-cavity mode, 
cii n is the input mode at the output port, and n is the 
leakage rate. This treatment of the transmission line res- 
onator as a cavity follows from the quantization of such 
quasi-one-dimensional transmission lines [Trl ] . The mode 
operator a of the cavity is related to the voltage and cur- 
rent carried by the transmission line by |l6j |: 

V(x,t) = ^[a(t)e ik * + a(t)U- ik *] 



I(x,t) 



V(x,t) 



(5) 



where L and C are the inductance and capacitance per 
unit length of the transmission line, and I is the length of 
the transmission line. For a given frequency, u>, the prop- 
agation constant k is fixed by k = cuV LC. These rela- 
tions come from solving the classical equations of motion 
for a lossless transmission line and quantizing the result 
[T(| . Note that once a(t) is given the solution is com- 
pletely specified, and therefore we characterize the state 
of the transmission line completely by defining a(t). 
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We will assume that a,i n (t) , the input mode at the out- 
put port is always the vacuum because the resonator can 
be engineered to minimize reflection and back-scatter. 
Given this model, we have the following equation of mo- 
tion for the intra-cavity mode |l5| : 



da(t) 
dt 



--[a(t),H ays ] - na(t) + \/~Kcn n {t) + \fKb in (t) 



-i[a(t), (u> r + — a z )a^a] — na{t) 



(t) + ^b in {t) 



(6) 



where H sys describes the dynamics of the intra-cavity 
mode, and bi n is the input mode at the input port of 
the resonator. Henceforth let x = 9 2 /^- If we Fourier 
transform this differential equation, use Eq. Q, and re- 
arrange, we can get a relationship between the input and 
output mode spectra: 



aWt(w) = — — ■ -(b m (uj) + a in (uj)) 

k + i(u r - lu + \a z ) 



- a, 



(7) 



where the a z in the denominator is simply a formal in- 
dication that the output spectra depends on the qubit 
state. So when the qubit is in a computational basis 
state [a z = ±1), the output mode at lu = uj r , the cavity 
resonance frequency, is: 



b in {u r ) T . a in (u> r ) (8) 



kHx 



Now, because a,i n (t) is the vacuum and the two input 
modes are assumed to be independednt, the second term 
of Eq. JSJ will not contribute to any normally ordered 
moments of a ou t, and we can effectively ignore it (even 
though it is strictly required to preserve commutation 
relations). Thus, 

2 

(i(Wr))Af = (—2— 2~frm(Wr))jV = ^ ^inM) (9) 

K + X 

where (-)jv' denotes a normally ordered moment, bi„ is 
a scaled version of 6j„ and ff^ is the qubit state de- 
pendent phase difference between the input and output 
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fields (tan(6' =F ) = T^)j exactly the quantity we want to 
measure in order to perform a QND measurement of the 
qubit. It can be measured using the standard heterodyne 
technique of mixing with a local oscillator of well defined 
phase to an intermediate frequency. However, in current 
experiments this mixing is done at room temperature af- 
ter amplification of the source signal, and this amplifi- 
cation and room temperature mixing both introduce too 
much noise for a high fidelity single shot measurement 
of the qubit state to be possible [TjI ■ It is precisely this 
problem that we address in this paper. In section llVl we 
propose two schemes for performing this QND measure- 
ment with high fidelity. 



III. MEASUREMENTS USING AN SET 

A single electron transistor (SET) is a device that is 
the combination of a small capacitance metallic island 
and high resistance source and drain leads. Electrons 
can tunnel from the leads onto the island and off again. 
Typically, one or more 'gates' are also capacitively cou- 
pled to the island, and voltages applied to these gates 
serve to change the electrostatic energy of the island and 
thus the electron tunneling rate. The typical SET op- 
erating regime is the so-called Coulomb blockade regime 
(which is set up by a choice of the source-drain voltage) , 
where the source-drain current is formed by single elec- 
trons tunneling between the contacts through the island 
which will only accommodate one excess electron at a 
time. In this regime, the voltage applied to the gate(s) 
has a dramatic effect on the tunneling rate and thus the 
device becomes a highly sensitive electrometer. 

In addition to the Coulomb-blockade regime, we can 
also operate an SET in the strong tunneling regime, 
where the source-drain voltage is larger than in the 
Coulomb-blockade regime, leading to larger tunnel con- 
ductances. In this regime the source-drain current is less 
sensitive to the gate voltage; in fact for near-optimally bi- 
ased source-drain voltage, the source-drain current versus 
gate voltage curve in this regime is well approximated by 
a sinusoid 0, ■ This is illustrated in Figure |3 We 
will sec that this strong tunneling regime is useful for 
performing measurements using an SET because of the 
well defined non-linear current-voltage behaviour. While 
operating in this regime, a DC bias voltage can be ap- 
plied to the gates to move the operating point along the 
voltage-current curve. The two most useful points along 
this curve from a measurement perspective are also illus- 
trated in Figure |21 



A. The SET as a mixer 

The non-linear relationship between the gate voltage 
and the source-drain current of an SET makes it possible 
to operate it as a frequency mixer, as demonstrated by 
Refs. [H|ll|l3. We shall use this feature of the SET 
extensively, and therefore in the following we will derive 
the exact relationship between the output current and 
gate voltage of an SET. 

Consider an SET operating in the strong tunneling 
regime with a voltage V g (t) = V g DC + Vi(t) + V 2 (t) ap- 
plied to the gate(s). We will analyze this situation be- 
cause it will be the relevant scenario when we consider 
making measurements in circuit QED using an SET (in 
that case V\ will correspond to the signal voltage from 
the cavity and V2 will be a local oscillator signal). Vq 
is a DC bias voltage that controls the SET operating 
point on the current-voltage curve of Fig. [5] We wish to 
obtain an expression for the source-drain current when 
the SET is operating at point A on the current-voltage 
curve. Because of the sinusoidal dependence on voltage, 



4 











Point A 






/ Foint B\ / 


\ / ^ 









FIG. 2: Approximate gate voltage versus source-drain current 
curve for an SET operating in the strong tunneling regime. 
The operating point of the SET is changed by applying a 
DC bias voltage to the gate. The two most useful operating 
points for measurement are denoted by Point A and Point B 
on the diagram. The x-axis is dimensionless and the y-axis is 
in arbitrary current units. 



we can write the source-drain current as: 
I(t) = Jo + A/qcos (^(V g DC + Vi(t) + V 2 (t))) (10) 



where V g is the DC bias voltage necessary to operate 
at point A. Now we will assume that V\ (t) and V 2 (t) each 
have a well-defined single frequency, and expand each out 
in a quadrature representation: 

I(t) = 7 + AJ cos[^(^ c + 

X 1 (t) cos(Qit) + Y 1 (t) sin(Qii) + 

X 2 (t) cos(Sl 2 t + Y 2 (t) sin(0 2 t + <£))](11) 

where 4> is the phase difference between V± and V 2 . 
Henceforth, let c g = 2irC g /e. Now we will use the sum 
angle formulae for cosine and sine to simplify. Firstly. 



I(t) = I + AJ [ cos[c g V g DC ] cos[c ff (Jfi(t) cos(Qif) + Yi(t) sin(fiii) + X 2 (t) cos(Sl 2 t + 0) + Y 2 (t) sm(Sl 2 t + 0))] 
- sm[c g V g DC ] sin[c fl (Jfi(t) cos(Qii) + Yi(t) sm(fM) + AT 2 (t) cos(ft 2 t + 0) + r 2 (<) sin(n 2 f + 0))]] 
= J + AJ /(i). (12) 

Now, because V g D is selected to be at point A, cos[c g V g DC ] = 1 and sm[c g V g ) = 0. This simplifies the fluctuation 
term in the above expression to: 

f(t) = cos[c 9 Xi(t) cos(Oit) + c g Yi(i) sin(f2it) + c g X 2 (t) cos(Sl 2 t + 4>) + c ff F 2 (t) sin(fi 2 i + 0))] 
= [ cos(c g Xi cos Slit) cos(c g li sin Slit) — sin(c fl Xi cos Slit) sin(c g Yi sin Slit)] ■ 

[ cos(c g X 2 cos(Sl 2 t + </>)) cos(c g y 2 sin(Sl 2 t + cf))) — sin(c ff X 2 cos{Sl 2 t + <j>)) sin(c g F 2 sin(fi 2 i + 0))] 
— [ sin(c g Xi cos Slit) cos(c g Yi sin Slit) + cos(c g Xi cos Slit) sin(c s Yi sin Slit)] ■ 

[sin(c 9 AT 2 cos(Sl 2 t + cf))) cos(c g Y 2 sin(Sl 2 t + (/))) + cos(c g X 2 cos(Sl 2 t + cf))) sm{c g Y 2 sin(f2 2 i + 0))] (13) 

where we have simply used the sum angle formulae repeatedly, and in the second line made the time dependence of 
the quadratures implicit. We want to analyze this expression in terms of its harmonics, and to do so we will expand 
the nested trigonometric functions in a Fourier-Bessel expansion pp| | using: 

oo 

cos(zcos(9) = J o (z) + 2^(-l) fc J 2fe (;z)cos(2fc60 

k=l 

oo 

cos(zsin#) = Jq(z) + 2^ J 2 k{z) cos(2fc6*) 
fe=i 



DC 



sin(zcos0) = 2^(-l) fc J 2fe+ i(z)cos((2£; + 1)0) 

fe=0 

oo 

sin(z sin 9) = 2 ^ J 2k +i{z) sin((2fc + 1)9) (14) 

fc=0 

where J„ is the n th order Bessel function. Note that when we apply this expansion to Eq. I|13J) the arguments to 
the Bessel function will be the quadratures of the two voltages; therefore we will get powers of Xi(t) and Yi(t) to all 
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orders. However, in applications of interest we will assume that these quadratures are small (i.e. C g V g <C e), and thus 
ignore contributions from higher powers. Since for small z, J n (z) — ► z n /2 n n\ |20| . we can truncate the expansions at 
J\ if we only consider first order contributions from the quadratures. Doing this yields: 

f(t) ps McgXjJoicgYjJoicgXdJofoYa) 

-[2Jx(c g X{)Jo(c g Yi) cos(Qii) + 2J (c g X{)J 1 (c g Yi) sin(fiii)] • 
[2J 1 (c ff X 2 )J o (c 9 Y 2 )cos(fi 2 t + 0) + 2J o (c 9 X 2 )J 1 (c 9 y 2 )sin(n 2 t + 0)] (15) 

Here we have also ignored the terms oscillating at 2f2i and 2f2 2 because in applications of interest these high frequency 
terms will not be detectable. Now expanding Jo and J\ to first order in their arguments, 

f(t) ps 1 - [2c 9 Xi cos(nii) + 2c g Y x sin(fiii)] • [2c 9 X 2 cos(fi 2 t + 0) + 2c 9 Y 2 sin(f7 2 i + 0)] (16) 

Finally, going back to the expression for the source-drain current, we get 

I(t) ps J + AJ [1 - 2 c 2x 1 X 2 (cos((0 1 -Q, 2 )t-4>) + cos((f2 1 +ft 2 )t + </»)) 
+ 2 C 2x 1 y 2 (sin((^ 1 - fi 2 )t - 0) - sin((Qi + fl 2 )t + 4>)) 

- 2c 2 g Y 1 X 2 (sm({{l 1 - n 2 )t -((>)+ sin((ni + fla)* + 0)) 

- 2 C 2y 1 r 2 (cos((ri 1 - n a )t - $ - cou((fii + n 2 )t + </>)) ] (17) 

Now if we ignore the high frequency terms (in applications of interest these frequencies will be above the bandwidth 
of the intrinsic SET response), we obtain the final expression for the source-drain current: 

I(t) ps I - 2A/ Cg [ X Y X 2 cos((f2i - Q 2 )t ~ <t>) ~ X X Y 2 sin((fii - Q 2 )t - </>) 

+Y X X 2 sin((f2i - Q 2 )t -</>) + Y X Y 2 cos((r2i - S7 2 )* - <P) ] (18) 

where Jo = Iq + Aio is the DC component, and the remaining signal is at the fl\ — f2 2 sideband frequency. This is the 
expression for the source-drain current at low frequencies and under the assumption that the gate voltage is small. 
Note that the SET has behaved as a mixer to downconvert the gate voltage signals to the intermediate frequency 
|f2i — 2 |. We will use this property in the next section when we examine high fidelity measurement strategies for 
circuit QED. 



B. SET measurement bandwidths 

As mentioned above, under the appropriate operating 
conditions the SET has a very high charge sensitivity, 
and this has led to several applications for it as an elec- 
trometer HUH miH. However, the SET has severe 
limitations when it comes to the measurement of high 
frequency signals. The intrinsic bandwidth of an SET is 
high as it is set by the tunneling rate of electrons thr ough 
the junctions - this is typically higher than 10 GHz |19l ]. 
However, the output bandwidth of an SET is constrained 
to the order of kilohertz because of the high resistance 
and the parasitic capacitance of the cabling that brings 
the SET output to room temperature [24[ . Typical values 
for the lead resistance and cable capacitance to the room 
temperature amplifier are 100kf2 and InF, thus giving an 
RC time constant of approximately 10 -4 , which is much 
too large to carry radio frequency signals without distor- 
tion. This bandwidth limitation of the SET makes it ill 
suited as a measuring device for circuit QED because of 
the high frequencies of the signals employed in the archi- 
tecture: the carrier frequency is the resonance frequency 
of the microwave cavity, about 10 GHz; and the band- 
width of the signal can be as much as the vacuum Rabi 



frequency of the qubit, about 100 MHz Q. Thus the 
output bandwidth of a typical SET is poorly matched to 
circuit QED frequencies. 

There are two well-known strategies to overcome the 
bandwidth limitations of an SET. The first is to use the 
non-linear current-voltage response of the SET to down- 
convert the signal from a large carrier frequency to an 
intermediate frequency by mixing with a local oscilla- 
tor |l7j . This allows measurements of signals at a wide 
range of carrier (center) frequencies, but the measure- 
ment bandwidth is still limited by the limited dynamic 
range of the output leads. 

The second strategy is to embed the SET in a high 
frequency resonant tank circuit, to create what has been 
dubbed the RF-SET, and measure the reflected (or trans- 
mitted) power from this tank circuit |25|. In this con- 
figuration, the SET acts as a variable resistor within an 
RLC circuit and the SET gate voltage effects the reflected 
power from the circuit by modifying the damping param- 
eters of the tank circuit. By designing the tank circuit 
to impedance match the SET and output lines, the SET 
measurement bandwidth can be increased by more than 
two orders of magnitude to over 100 MHz |25|]. Due to 
the mixing properties of the SET, the spectrum of the 
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reflected power from the tank circuit has components at 
the intermediate frequencies \f g ± /lcI where f g is the 
frequency of the gate voltage and Jlc is the resonance 
frequency of the tank circuit. Thus the RF-SET con- 
figuration has the dual effect of increasing the output 
bandwidth, and downconverting a high carrier frequency 
signal to an intermediate carrier frequency. But note that 
the downconvcrsion in this case is to a fixed intermediate 
frequency (set by the tank circuit parameters) unlike in 
the previous case where the intermediate frequency could 
be changed by the local oscillator frequency choice. 

Recently, Swenson et. al have combined both of these 
innovations, the local oscillator mixing and the RF-SET, 
to create an electrometer that has a large output band- 
width as well as a large, tunable center frequency [T3 |. In 
the next section we will describe two ways in which the 
electrometer of Swenson et. al can be used to perform 
measurements in circuit QED. 



IV. TWO PROPOSALS FOR HIGH FIDELITY 
MEASUREMENT 

In this section we propose two strategies for performing 
high fidelity QND measurements of the qubits in circuit 
QED. Both schemes involve using an RF-SET as a solid- 
state mixer that can be fabricated onto the same circuit 
as the qubits and waveguide. From Sec. |II|we see that 
the microwave field transmitted from the cavity contains 
information about the state of the qubits because of the 
intra-cavity coupling. However, using this transmission 
line directly as a gate to an RF-SET would not give us the 
information we need because that would simply measure 
the intensity of the field, and the qubit state is encoded 
in the amplitude and phase of the field. This is a stan- 
dard problem in quantum optics (because photodctectors 
detect intensity) and there the solution is to mix the sig- 
nal with a local oscillator. The same technique is used 
in circuit QED 0, however the mixing is done at room 
temperature. In the following we describe how to do this 
mixing on the circuit using an RF-SET. 



A. Scheme 1: Using a dual-gate SET 

We can perform the mixing with a local oscillator us- 
ing a dual gate SET at the core of the RF-SET. Figure 
13 is a schematic of a dual gate single electron transistor. 
We will assume the source-drain voltage is set such that 
the SET is in the strong tunneling regime. To perform 
heterodyne mixing using this device we apply the signal 
voltage, V s of frequency uj r , to one gate and a local oscil- 
lator, Vlo of frequency ujlo, to the other gate. Then the 
source-drain current is modulated by both gate voltages. 
We also bias the SET such that it is operating at point 
A of Fig. [5] Then the source-drain current is: 

I(t) = I + AI cos (^(V g DC + V S + V LO )) (19) 



Then under the assumption that both V s and Vlo are 
small, we can follow the calculations in Sec. Illll to obtain 
an expression for this current at the sideband frequency 
uo r — luloi m terms of the quadratures of the voltages: 



I(t) 



1tt 2 AI C% r 

5 [X s Xlo cos((w r - w L o)t 



- X s Y L o sin((w r - Lo L o)t 
+ Y s X L o sin((w r - LULo)t - 
+ Y s Y L o cos((w r - uJ L o)t - 



(20) 



where X s , Y s (Xlo, Ylo) are the quadratures of the sig- 
nal (local oscillator), and <j> is the phase difference be- 
tween the two voltages. In addition, because we are 
using an RF-SET the output signal has components at 
|aJic±(a; r — (jJlo)\ where ujlc is the tank circuit resonant 
frequency. At this point we want to choose lolo ~ w r and 
thus have an output signal around the ujlc center fre- 
quency. Therefore let lolo = w r and consider the phase 
choice <j> = 0. Then considering just the change in the 
current from the static value (Iq) we have 



AI(t) 



8ir 2 AI C 2 



[ X s X LO + Y s Y LO ] (21) 



Expanding the quadrature operators in terms of mode 
annihilation and creation operators (X s = a + a\Y s = 
i(a — a'), Xlo = b + b^ , Ylo = i{b — w), where a is the 
cavity output mode annihilation operator whose explicit 
expression in terms of the input mode was derived in 
section [n] - i.e. a ou t in Eq. ©) and considering the 
expectation value of this current under an arbitrary state, 
we have 



(AJ(t)) 



167T 2 A/nC2 



?-{a)b + tfa) 



(22) 



at bj — ujlc- Now finally, if we prepare the local oscilla- 
tor in a coherent state \/3) (for j3 real because the relative 
phase difference between the local oscillator and the sig- 
nal has already been accounted for by the choice of <f>) 
then we have 



(AI(i)> 



Wn 2 AI C 2 g 



,t\ 



(23) 



That is, the change in the average monitored current is 
proportional to the X quadrature of the signal. Similarly, 
if we had chosen the phase to be (j> = tt/2 the average cur- 
rent would have been proportional to the Y quadrature 
of the signal voltage. 

Using such a mixer with the source signal being the 
transmitted field from the circuit QED transmission line 
will allow one to detect the phase shift of the field mode, 
and thus the state of the qubits. The advantage of such 
a scheme is that measurement apparatus, in this case 
the RF-SET and the local oscillator, can be fabricated 
directly onto the same circuit as the qubits and resonator 
and operated at the cryogenic level (e.g. Hi|). This 
would lead to a higher fidelity measurement signal. 
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Gate I 



Source Drain 




Gate 1 



FIG. 3: Schematic of a dual gate SET. The black box in 
the middle is the island which is capacitively coupled to the 
source and drain leads and also the gates. The gate voltage 
for a dual gate SET is simply the sum of the voltages applied 
to the two gates. 




put signal of the four-port hybrid is made the gate signal 
of a single gate RF-SET operating in the strong tunnel- 
ing regime and at point A on Figure |2 (i.e. the DC bias 
voltage is selected so that the current- voltage curve is at 
an extremum). This situation then becomes exactly the 
same as a dual gate SET with the source and local oscil- 
lator each being applied to separate gates. In this case V\ 
takes the place of the signal voltage, V2 takes the place 
of the local oscillator, and as before the sum of the volt- 
ages modulates the source-drain current. Therefore this 
setup becomes entirely analogous to the RF-SET mixer 
measurement of section TlVAl and the analysis from that 
section can be repeated: if we choose the frequencies of 
Vi and V2 to be equal, and both voltages to be small, 
quadrature measurements of the cavity output mode can 
be made at the tank circuit resonant frequency u>lc- 

The fabrication of SETs with multiple, well character- 
ized gates is a challenging exercise, and therefore this 
scheme will be easier to implement than the scheme 
of section II V Al if the hybrid coupler is easier to engi- 
neer than a dual-gate SET. Issues such as cross-coupling, 
losses and asymmetries at the SET gates are less impor- 
tant with this second scheme. The ideal hybrid coupler is 
a lossless, unitary device that is likely to be much easier 
to characterize than a multi-gate SET. 

Finally, we note that Mariantoni et. al. have recently 
proposed the use of an on-chip hybrid ring coupler to per- 
form homodyne measurements in a manner very similar 
to the method above [30l| . 



FIG. 4: Schematic of the transmission line coupling that mim- 
ics the behaviour of a beam-spiltter in optics. Vi and V2 are 
the input signals and the outputs are the sum and difference 
of these. 

B. Scheme 2: Engineering a mixer from waveguides 

Our second scheme for achieving a high fidelity QND 
measurement of the qubits involves engineering a copla- 
nar waveguide coupling to mimic a quantum optics beam- 
splitter. A schematic of such a setup is shown in Figure 
0J The box in the middle generates the sum and differ- 
ence of the signals provided as input. In the circuit QED 
model, all the signals are carried on coplanar waveguides 
and therefore the thick lines represent these waveguides 
while the arrows indicate the direction of signal propaga- 
tion. Such a coupling can be engineered through the use 
of well known four-port hybrid devices, for example, the 
rat-race hybrid |27ll28ll23 |. One of the input waveguides 
will be the transmission line that couples to the qubits 
(i.e. the 'cavity'), and the other will be a new coplanar 
waveguide that carries a local oscillator mode. 

Given such a coupling, we can use a single-gate RF- 
SET to perform intensity measurements of the sum (or 
difference) signal. That is, the sum (or difference) out- 



V. A FEEDBACK PROTOCOL 

The availability of high fidelity measurement enables 
the implementation of quantum feedback protocols, and 
in this section we describe one simple feedback protocol 
to deterministically generate and stabilize entanglement 
between two qubits in the circuit-QED architecture. 

As mentioned in section [n] in the dispersive regime 
the interaction between the microwave cavity mode and 
a qubit takes the form Hi nt = (hg 2 / A)a z a ji a. If we have 
two qubits both coupled to the same microwave mode 
through the Jaynes-Cummings interaction, then the ef- 
fective Hamiltoninan in the dispersive regime is well ap- 
proximated by 0, : 

Hdisp « (uj r + 2xJ z )a^a+ (uj a +x)Jz 

+xKcr 2 -+crr4) ( 24 ) 

where h has been set to 1, J z = \(a\ + erf), x = 9 2 /A> 
and af and a~ are the qubit raising and lowering oper- 
ators for qubit i. The last term is an exchange coupling 
between the qubits that is induced by their interaction 
with a common mode. It is convenient to work in an 
interaction picture and simplify the Hamiltonian to: 

Hint ~ 2xJ z a^a + x(<?t a 2 + a \ a i) ( 25 ) 
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The qubit-ficld interaction can be viewed as one where 
the field measures a collective property of the qubits. It 
has been noted by several authors that such an interac- 
tion combined with a monitoring of the output field can 
probabilistically generate an entangled state between the 
two qubits |3^ . Isifl . l34l | . Additionally, in a series of elegant 
papers, Mabuchi et. al. have shown that it is possible to 
deterministically generate entanglement in a cavity QED 
syste m ( with a similar Hamiltonian) by employing feed- 
back [35l l36l I37J . We will show that it possible to do the 
same in circuit QED. 

The motivation for introducing feedback in order to 
generate and stabilize entanglement comes from the fol- 
lowing simple observation. Consider the evolution of an 
initial product state of two qubits polarized in the x di- 
rection and a coherent state of the field: 

(|00) + [01) + 1 10) + |11»® |a) 
-> |00) \ae~ l2xt ) +e- lx *|01) \a) + 

e- lxt |10)|a) + |ll)|ae l2xt ) (26) 

where the arrow indicates unitary evolution for time t 
under the interaction Hamiltonian Eq. i|25|) . Notice that 
the field acquires a phase shift if the qubits are in the 
|00) or 1 11) state. Thus, if we apply a J x = \{o\ + erf) 
rotation when we see a phase shift in the field, it is pos- 
sible that the system could be forced into the entangled 
state \<f> + ) = ^(|01) + 1 10) ) . We will show that this is 
indeed the case. 

The first step in modeling the dynamical system is to 
include the driving and damping of the microwave cavity 
mode. This results in the following master equation for 
the system: 

^ = -i[H,p] + KV[a]p (27) 

where k is the rate of damping of the cavity, and the 
Hamiltonian with the added driving term is 

H = e(a + a f ) + x(2J 2 a f a + af + 0^0^) (28) 

Now if we perform a homodyne measurement of the cav- 
ity field in order to detect the phase shift, then in analogy 
with the one qubit case, we gain information about the 
J z eigenvalue of the two-qubit state. The measurement 
current increment is expressed as 0, 0, l3Sj 

dl(t) = n(a + a>)dt + ^KdW{t) (29) 

where dW(t) a random variable that represents the 
noise; it is a Wiener increment: E[<W(t)] = 0, and 
E[dW(i)<iW(s)] = S(t—s)dt 1 . If we assume the measure- 
ment is unit efficiency, we can represent the evolution of 



the system as a pure state evolution (also called an unrav- 
elling of Eq. (I27[l into a stochastic Schrodinger equation 

(sse)) a a® 

d \tp(t)) = [-W - |a f a] \ij)(t)) dt + dl(t)a \tp(t)) (30) 

This is an unnormalized version of the SSE, and H the 
Hamiltonian is from Eq. I|28|l . This is a filtering equation 
that modifies the state evolution to include the back- 
action of the measurements while also enforcing consis- 
tency with the observed measurement results. 

Now we need to add the feedback. The simplest form 
of feedback is direct Hamiltonian feedback [3!j, where 
the feedback operator is conditioned by a signal that 
is linearly related to the current (and only the current) 
measurement value. Such feedback results in Markovian 
system dynamics and thus an evolution equation that is 
particularly simple. However, this type of feedback is 
poorly suited to tasks such as stabilization because of 
the amount of noise in the instantaneous measurement 
signal. Therefore, in order to deterministically gener- 
ate and stabilize entanglement in our system we will use 
a state-estimate based feedback mechanism that results 
in non-Markovian, but generally more stable system dy- 
namics. Essentially the task is to estimate the J z value of 
the qubit state with hig h accuracy and base the feedback 
upon it. In Refs. [3|| |32j Mabuchi et. al. investigate 
the optimality and robustness of such an approach in 
achieving the same goal in a similar system and conclude 
favorably on it. We will not do the optimality analy- 
sis, but will simply illustrate the validity of the feedback 
protocol by simulation. 

We already have a noisy, continuous measurement of 
J z (the homodyne current measures a phase shift, which 
is directly proportional to the J z value of the qubits), and 
we will refine it into a less noisy estimate of (J z ) by two 
steps: (1) low-pass filter the measurement signal over a 
small time window, and (2) condition the feedback with a 
power of the filtered measurement signal. The first step 
forms a crude, but efficient estimate of the observable 
(Jz), and the second reduces the noise in this estimate 
before conditioning the feedback with it. In the first step 
the instantaneous measurement is filtered using a simple 
low-pass filter to produce the smoothed signal 

m = h f t T e " 1{t ' t ' )dI[t ' ) (31) 

where the factor Af normalises the smoothed signal to an 
maximum magnitude of unity. We model the feedback by 
simply adding a conditioned Hamiltonian evolution term 
to the SSE of Eq. ® 

d\ip(t)) = [-iH - ^a}\ip(t)}dt + dl(t)a\ip(t)) 

- i\R(t) p J x \^(t)) dt (32) 



1 E[-] denotes the classical expectation value. 



where P is the power to which the smoothed sig- 
nal is raised before conditioning the feedback, and A 
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parametrises the feedback Hamiltonian strength. Such 
feedback could be performed in the circuit QED architec- 
ture by utilizing the electronic control over the Cooper 
pair box qubits, or by irradiating the cavity with an ad- 
ditional mode in resonance with the qubit transition fre- 
quency as done in Ref. [l2^ . 

At this point, we note that there arc a number of free 
parameters in this SSE: e, the field driving strength; \i 
the qubit-field coupling strength; k, the field damping 
rate; A, the feedback rate; T and 7, the parameters of 
the low pass filter; and finally, P, the power to which 
the smoothed current is raised. We will illustrate the 
effectiveness of the above feedback protocol in generating 
a maximally entangled two qubit state by numerically 
solving the (non-linear, non-Markovian) SSE, Eq. I|32|) . 
for the following parameter values: e = 100, \ = 25, k = 
100, A = 100, P = 3. The time step is dt = 0.0001, the 
filter values are T = 2000 x dt and 7 = 0.003, and the 
total time simulated is 10 time units (the units of time 
are arbitrary) . The choice of these parameters were fairly 
arbitrary, however, there are some guiding principles for 
the relationship between the parameters for the scheme 
to be effective: 

1. We want a heavily damped cavity so that we have 
as much information about the qubits as possible. 
We also want to minimize the long-term entangle- 
ment between the qubits and the cavity mode be- 
cause such entanglement would decrease the entan- 
glement between the qubits themselves. Essentially 
we want the mode to act like an ideal measurement 
device. These consideration suggest e k k » ^. 

2. The feedback should be strong enough to compete 
with the damping effects of the field. Thus A > k. 

3. The filter should be large enough to reduce the 
noise to acceptable levels but still be responsive to 
signal changes. 

4. We should raise the smoothed signal to a power 
such that the noise is reduced while the signal is 
still preserved. The quantity P = 3 was arrived at 
by trial. 

Also, we truncate the harmonic oscillator state space to 
25 dimensions for the numerics. This truncation is valid 
because the strong damping keeps the oscillator state 
close to the vacuum. The stability of the results derived 
using this truncation were confirmed by comparing to 
results using a larger state space. 

Solving the SSE in Eq. (|32|l once is not sufficient to 
make conclusions about the properties of the feedback 
protocol because that would only determine the perfor- 
mance under one noise record {dW(t)}. We are more 
interested in the average properties of the entanglement 
generation scheme, and thus the SSE is solved 300 times - 
300 trajectories are generated - all with the initial state 
\i>o) = (1(100) + |01) + 1 10) + |11)) ® \a) with a = 3. 
Then the average entanglement generated is used as a 
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FIG. 5: The evolution of entanglement averaged over 300 
trajectories. The parameters used in the simulation are in 
the main text. The time units are arbitrary. 
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0.4 V 
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time 

FIG. 6: The evolution of fidelity with the maximally entan- 
gled state \(f> + ) = ^TjdOl) + 1 10)) averaged over 300 trajecto- 
ries. The parameters used in the simulation are in the main 
text. The time units are arbitrary. 



measure of performance. Figures |S] and [S] summarize 
the results of these simulations. Figure |3] shows the 
evolution of the average entanglement (measured using 
the concurrence measure |4*cj |) of the two qubits, and 
figure HJ| shows the evolution of the average overlap of 
the two qubit state with the maximally entangled state 
\<j) + ) = ^(|01) + 1 10)). These figures clearly indicate 
a convergence to the desired entangled state, in fact we 
find that 99% of the 300 trajectories converge onto the 
\<fr + ) state. Furthermore, as t — + 00, the convergence 
percentage approches 100%. It is important to note that 
the measurement record indicates the state of the qubits, 
and therefore we can tell when we have converged on the 
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entangled state. 

Thus we have convincing numerical evidence to sug- 
gest that the above feedback scheme will converge on the 
maximally entangled state \4> + ) in a short time with very 
high probability. 

VI. CONCLUSION 

We have presented two schemes for achieving high fi- 
delity single-shot measurements in the circuit QED archi- 
tecture. As with original proposal Q, the measurement 
is done in the Josephson basis and therefore minimizes 
excursions of the CPB from the optimal charge degen- 
eracy point. Furthermore, the use of dispersive readout 
and an RF-SET minimizes on-chip power dissipation, an 
important consideration for solid-state implementations 
of quantum computing. Both measurement schemes use 
the RF-SET as a fundamental element, and the possi- 
bility of making them quantum-limited hinges crucially 
on the performance of these SETs. It is conceivable that 
the ultimate limit to the performance of these measure- 
ment schemes is determined by the shot-noise through 
the SETs. We briefly note that it may be possible to 
use a Josephson junction as an alternative to the SET 
to perform the measurements. The sinusoidal current- 
flux relationship of a large capacitance Josephson junc- 
tion (Ej > E c ) [13 could be used to generate the non- 
linearity needed to mix the signal and local oscillator 
voltages. 

Despite the similarity in the analysis of the two mea- 
surement schemes, they are very different physically. 



Both bear a resemblance to quantum optics heterodyne 
measurement, with the second being almost the solid- 
state counterpart of it. However, a major departure from 
the quantum optics analogy is the fact that for these SET 
mixer strategies, we need the signal and the local oscil- 
lator amplitude to be small. This is a consequence of the 
ability of SETs to generate non-linearities of all orders. 

High fidelity measurements in the circuit QED archi- 
tecture will enable many interesting experiments in quan- 
tum feedback control, and we have briefly discussed one 
such experiment to generate and stabilize entan glem ent 
between two qubits. The results of Refs. [35l l36t 13^ 
suggest that multi-qubit entanglement can also be gen- 
erated in circuit QED using a similar protocol to the one 
we have described. In addition to quantum control appli- 
cations, high fidelity single shot measurements arc likely 
to be crucial for quantum error correction, an essential 
component of quantum computing. 
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